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ABSTRACT 

In this paper we constrain the evolution in the comoving space density for low- 
luminosity (typically FRI) radio sources as a function of cosmic epoch by matching 
deep radio surveys (WENSS, FIRST and NVSS) with the Sloan Digital Sky Sur- 
vey (SDSS). This results in 1014 matched radio sources covering an effective area of 
0.217 sr, 686 with L325MHZ < 10 25 W Hz -1 sr _1 , which is an order of magnitude larger 
than any previous study at these luminosities. Using the non-parametric V/V m:ix test 
we show that low-luminosity radio sources evolve differently to their more powerful, 
predominantly FRII, counterparts. Indeed, we find that the lower luminosity popula- 
tion is consistent with a constant comoving space density with redshift, as opposed to 
the strong positive evolution for the more powerful sources. 
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1 INTRODUCTION 

The evolution of radio-loud active galactic nuclei (AGN) has 
been a difficult quantity to pin down for all but the most 
powerful sources up to moderate redshifts. These sources ap- 
pear to trace the most massive galaxies at all cosmic epochs 
. (e.g. Jarvis et al., 2001; Willott et al., 2003) and so under- 
standing their evolution is of strategic importance to astron- 
omy. Further, radio sources make particularly useful probes 
of galaxy and AGN evolution because radio emission is not 
affected by dust obscuration. The evolution in the space den- 
sity of the radio source population is therefore a crucial el- 
ement in our understanding of the AGN phenomena. 

It has been suggested that the radio AGN population 
is composed of different types of objects. The more pow- 
erful sources show a high surface brightness at the ends 
of a double lobed structure (Fanaroff-Riley type II; here- 
after FRII; Fanaroff & Riley 1974) whereas the less power- 
ful sources have high surface brightness nearest their cen- 
tres (FRI). It is now well established that the FRII radio 
sources evolve strongly between z ~ 2 — > (Longair 1966; 
Dunlop & Peacock 1990) with a decrease in the comoving 
space density of ~ 1000 since z ~ 2. However, very lit- 
tle is known about the lower- luminosity FRI radio sources. 
Various groups have tried to understand how the lower- 
luminosity radio sources evolve with redshift. The conclu- 
sions drawn from these studies vary. For instance, Jackson 
& Wall (1999) find evidence for possible negative evolution. 
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Laing, Riley & Longair (1983) find, from 17 sources in the 
complete 3CRR sample, evidence for no-evolution. Similarly, 
a study by Waddington et al. (2001) of 72 sources also sug- 
gests little evidence of evolution in the space density between 
z = and z = 1 in the low-luminosity regime. However, 
Snellen & Best (2001) find, from two objects in the Hub- 
ble Deep Field and the flanking fields, evidence for positive 
evolution. This fundamental discrepancy arises because of 
the observational difficulty in establishing a large complete 
flux-density-limited radio sample at faint fluxes due to the 
vast overall increase in radio source counts. 

The problem of small number statistics can now be over- 
come with large multiwavelength surveys. The advent of the 
large radio surveys such as the 325 MHz Westerbork North- 
ern Sky Survey (WENSS); the Faint Images of the Radio 
Sky at Twenty-cm (FIRST) and the 1.4 GHz NRAO VLA 
Sky Survey (NVSS) along with large scale optical surveys 
such as the Sloan Digital Sky Survey Data Release 1 (SDSS 
DR1) allow us to probe the low-luminosity radio source pop- 
ulation in a statistically complete manner. These surveys, 
when combined, provide the most efficient method of inves- 
tigating the evolution of the faint radio source population as 
a function of cosmic epoch. The limiting factors are the flux- 
density limits of the radio surveys and the depth of the opti- 
cal data. The WENSS and FIRST surveys can probe a radio 
source with L325MHZ = 10 25 W Hz -1 sr -1 (roughly at the 
FRI/FRII break luminosity) at a flux limit of 5*325 = 40 mjy 
up to z ~ 1. The most limiting constraint comes from an 
upper limit in the optical flux from the SDSS (we adopt 
i < 21). 

In this paper we investigate the evolution in the co- 
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moving space density of low- luminosity radio sources with 
cross-matched catalogues formed from the WENSS, FIRST 
and NVSS radio surveys and the recently released photomet- 
ric survey SDSS DR1. In §2 we describe the selection of the 
radio sources and their matching to the SDSS DR1. §3 con- 
tains details of the calculation of the photometric redshifts. 
In §4 the evolution in the comoving space density of radio 
sources is tested via the V/Vmax statistic. In §5 we present 
the results and §6 is a summary of our conclusions. 

We define radio spectral index a r ad as S v oc j/ -a!rad , 
where £>„ is the flux-density at frequency v. We assume 
throughout that H B = 70 km s~ 1 Mpc^ 1 and a Qm = 0.3, 
Aa =0.7 cosmology. 



2 SAMPLE SELECTION 

We select a sample area that is sufficiently large to in- 
crease the number of low-luminosity, predominantly FRI, 
radio sources from ~ tens to about 1000. We choose the 
continuous area of the northern SDSS DR1 sample which is 
bounded by FIRST at a declination of 64° and WENSS at 
38°. The shaded are in Fig.^shows the resulting area of our 
survey, covering 713 square degrees (0.217 sr). 

2.1 Radio source selection 

We wish to select radio galaxies where the optical light is 
dominated by the host galaxy, rather than quasars, and the 
radio emission may undergo Doppler beaming effects along 
the line-of-sight (e.g. Jarvis & Rawlings 2000). Therefore, 
we start with the low-frequency WENSS radio survey that 
will preferentially select optically thin steep spectrum lobe 
emission rather than optically thick flat spectrum (possibly 
beamed) emission. 

The WENSS radio sources are matched to NVSS and 
FIRST radio surveys, bounded by the SDSS DR1 photomet- 
ric survey. We summarise each of these surveys only briefly 
here - details can be found at each of their respective web- 
sites and references therein. 

The WENSS survey (Rengeli nk et al. 1997; 
http:/ /www. strw.leidenuniv.nl/wenss/ I covers 3.14 sr north 
of +30° declination at 325 MHz with resolution of 54" and 
a flux-density completeness limit of approximately 18 mjy 
(5<t). We choose a conservative flux-density limit of 40 mjy 
which provides the limit for this survey. The NVSS radio 
survey (Condon et al. 1998; http://www.cv.nrao.edu/nvss/ ) 
is a 1.4 GHz continuum survey covering 10.3 sr of sky north 
of -40° declination. The completeness limit is about 2.5 mjy 
and the images have a 45" resolution and a nearly uniform 
sensitivity. The FIRST survey (Becker, White & Helfand 
1995; White et al. 1997; http://sundog.stsci.edul is at 1.4 
GHz and covers 2.75 sr at -10° < declination < 64° with 
resolution of 5" and a 1 mjy source detection threshold. 

2.2 Radio source matching 

We cross-match the WENSS and NVSS catalogue positions 
taking all matches with positional offsets of ^ 40" as can- 
didate radio detections. (We use the low-resolution NVSS 
survey over FIRST to determine 1.4 GHz flux-densities as 
some of the flux may be resolved out by the higher-resolution 
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Figure 1. The survey area of the matched radio and photometric 
sample covering a total area of survey is 713 square degrees (0.217 
sr). The curves are the area covered by the SDSS DR1 which are 
bounded by FIRST at a declination of 64° and WENSS at 38° . 
There are 1014 matched a ra( j > 0.5 radio sources in the survey 
(shown as dots). 

FIRST survey.) Increasing the distance to ^ 45" increases 
the number of matched targets by only ~ 1%. At this dis- 
tance the probability of choosing a NVSS source within an 
arbitrary position is 2% whereas for the WENSS catalogue 
this is 0.8%, so we expect the contamination to be at the 
level of ~ 1% and the 40" offset is considered optimal. 

We calculate the radio spectral index, Q ra( j, for each 
source using the flux-density at frequencies 325 MHz and 1.4 
GHz from the matched radio sources in the low-resolution 
WENSS and NVSS surveys respectively. In order to select 
against quasars and other non-lobe dominated sources we 
choose a ra d > 0.5. 

The NVSS coordinates of this cross-matched sample are 
further matched to the high resolution radio galaxy cata- 
logue, FIRST. Again we choose a positional offset of ^ 40". 
We do this for two reasons: (i) the increased resolution of 
FIRST allow us to improve the positional accuracy of the 
matched WENSS-NVSS catalogue, thus reducing the num- 
ber of mis-identified sources with the photometric survey 
and (ii) the FIRST survey image cutouts allow an accu- 
rate visual check on the matching procedure to the SDSS 
DR1. The matched WENSS-NVSS-FIRST yield a total of 
7898 matched radio targets in the sample area to be then 
matched to the photometric survey. 

2.3 Photometric matching 

We cross-match the SDSS DR1 catalogue with the FIRST 
matched positions, with positional offsets of ^ 5". We use 
the SDSS modified Petrosian magnitudes (Petrosian, 1976; 
Strauss et al., 2002). In this magnitude system the aper- 
ture size depends on the shape of the galaxy's radial surface 
brightness profile but not its amplitude. The disadvantage 
of this system is that the photometric measurements be- 
come increasingly noisy for r fainter than around 20. The 
other option is to use SDSS model magnitudes which are 
less noisy for the faint galaxies but have serious systematics 
for galaxies r < 20. The SDSS DR1 quote point spread func- 
tion magnitude limits (95% detection repeatability for point 
sources) of 22.0, 22.2, 22.2, 21.3, 20.5 in u,g,r,i,z respec- 
tively with the errors in the magnitudes ^ 3%. For Petrosian 
magnitudes the reported errors are considerably larger than 
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Figure 2. Photometric redshifts versus spectroscopic redshifts for 
265 radio sources in our sample with spectroscopic redshifts. The 
ANNz code is trained on ~ 120, 000 galaxies from the SDSS. A 
linear least-squares fit to the sample with z 0.5 gives Z p ^ ot = 
Z Bp cc * 0.943 — 0.003, shown as a solid line. The dotted line is 

■^phot — -^spec- 



this. To ensure our sample is complete, and to minimise the 
photometric errors, we take a colour cut of i < 21. This 
effectively limits our analysis to z ^ 1, assuming the host 
galaxies are typically ~ 3L* (McLure et al. 2004) - which 
makes it the most limiting factor in the survey. Additionally, 
we use the SDSS morphological classification and select only 
non-stellar sources thus selecting against quasars. 

We overplot the optical matches from SDSS DR1 onto 
the radio images from FIRST. In most cases the radio source 
is unambiguous. However, careful visual inspection of the 
images reveals that 8% of the sources with optical identifi- 
cations within 1' are not easily associated with a core of a 
radio galaxy. The majority of these cases are due to mul- 
tiple identifications along the radio axis, which presumably 
arises from the well known alignment effect of the optical 
and radio emission (e.g. Eales et. al. 1997; Best, Rottgering 
& Longair, 2000; Inskip et al. 2002). As the narrow-line lu- 
minosity is correlated with radio luminosity (e.g. Willott et 
al. 1999) then this effect will be most prominent in the more 
powerful FRII subset. Thus we are confident that the sur- 
vey as a whole is > 90% complete, and the low-luminosity 
sources - which are our main concern in this paper - have 
a higher completeness. These selections result in 1014 radio 
sources with associated photometric data. 

We incrementally alter the search radius of the optically 
matched FIRST and SDSS DR1 radio sources from 5 arcsec 
to 2 arcsec in order to investigate the effect on the contami- 
nation of mis-identified sources and the completeness of the 
sample. Ivezic et al. (2002) discuss this problem in detail. 
These authors find that with a 1 arcsec search radius the 
sample is 72% complete and 1.5% contaminated. Increasing 
the search radius to 3 arcsec results in > 99% completeness 
and 9% contamination. We expect the contamination in our 
survey to be significantly less than that found by Ivezic et 
al. as the matching of each radio source is visually inspected 
by overplotting the SDSS DR1 photometric data on maps 
of the FIRST radio sources, whereas the Ivezic et al. study 
was automated. 



3 PHOTOMETRIC REDSHIFTS 

Redshift information is crucial in determining the evolution 
in the comoving space density of radio sources. Gaining spec- 
troscopic redshift information for all of the radio sources in 
our sample would be extremely time consuming on any tele- 
scope. However, the five-band photometry from the SDSS 
DR1 allows us to use photometric redshifts to make rea- 
sonably accurate estimates of the whole sample. We use two 
publicly available codes to do this. The first is Hyperz (Bol- 
zonella, Miralles & Pell, 2000) , which uses galaxy templates 
and various degrees of reddening to determine the most 
likely redshift. We assume a priori that the host galaxies 
of our radio sources are either elliptical galaxies or bulge 
dominated spirals (e.g. McLure et al. 2004), and as such 
only use the corresponding template spectra. 

The second code (ANNz; Collister & Lahav 2004) esti- 
mates photometric redshifts by employing neural networks 
which require a large training set of galaxies for which the 
redshifts are already known. Collister & Lahav demonstrate 
the effectiveness of their method using SDSS DR1 and com- 
pare it to the Hyperz technique. They find that when a large 
training set is available the ANNz technique is more robust 
than Hyperz method with smaller random and systematic 
errors. 

We use the large sample of galaxies with spectroscopic 
redshifts and five-band photometry from the SDSS DR1, 
which comprises ~ 120, 000 galaxies. In Fig. H we show the 
accuracy of the ANNz technique for targets in our sample 
with spectroscopic measurements. The random error on the 
photometric redshifts is typically Az ~ 0.03 at z < 0.5. 
Inspection of Fig. [5] reveals a slight systematic error. We 
discuss the effect of this systematic underestimate of the 
photometric redshifts in f|H 

Beyond z > 0.5 the errors in the ANNz analysis become 
larger, predominantly due to the fact that the majority of 
the galaxies observed in the spectroscopic training sample 
have i < 19.5. The relative dearth of high-redshift sources in 
the spectroscopic training set fainter than this limit means 
that ANNz is not the ideal photometric redshift code for 
this region of parameter space. Therefore we use Hyperz 
estimates for sources with i > 19.5. 



4 THE V/Vmax STATISTIC 

In order to assess the evolution in the comoving space den- 
sity of radio sources we use the non-parametric V/V ma x 
method (Schmidt 1968; Rowan- Robinson 1968). This tech- 
nique is useful as we do not have to make any a priori as- 
sumptions as to the form of any evolution. It also allows 
us to easily incorporate additional selection criteria. Indeed 
this is critical for this study as we have an optical magnitude 
limit of i < 21.0 mag. 

Under the null hypothesis of a uniform distribution, 
the value of V/Vmax will be uniformly distributed between 
and 1. For such a sample the mean value is (V/Vmax ) = 
0.5±(12A r )~ 2 , where N is the number of objects in the sam- 
ple. A value of {V/Vmax ) > 0.5 indicates that the sources 
are biased towards larger distances, or an increasing space 
density with redshift; whereas {V/Vmax ) < 0.5 indicates a 
deficit of high-redshift sources, or a decline in the space den- 
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Table 1. 


V/Vmax 


for 8 bins 


in 325 MHz radio luminosity. 


^325MHz 


N 


V/V ma x 


la 


< z > 


V/Vmax (z ad j) 


23-23.5 


55 


0.43 


0.04 


0.11 


0.42 


23.5-24 


125 


0.42 


0.03 


0.18 


0.43 


24-24.5 


241 


0.52 


0.02 


0.31 


0.52 


24.5-25 


265 


0.45 


0.02 


0.45 


0.46 


25-25.5 


192 


0.54 


0.02 


0.66 


0.53 


25.5-26 


86 


0.62 


0.03 


0.79 


0.63 


26-26.5 


31 


0.70 


0.05 


0.96 


0.72 


26.5-27 


13 


0.62 


0.08 


0.91 


0.64 



sity evolution. The total fraction of sources for which Vmax 
is set by the optical selection criteria is 35%. This is split 
into 13% of the low-luminosity sources and 83% of the high- 
luminosity sources. To incorporate the additional optical se- 
lection into the V/V m ax statistic we need to know the redshift 
at which the galaxy fades to i > 21 mag. This is found by 
fitting the multi-colour data with various galaxy templates 
generated from the GISSEL96 synthetic galaxy spectra li- 
brary of Bruzual & Chariot (2003). We use ten different 
templates of an evolving elliptical galaxy, with the redshift 
given by the photometric redshift techniques described in 
§01 The spectral fit to the data is minimised with a simple 
X 2 fitting routine and subsequently inspected by eye (this 
also acted as a further consistency check of the photometric 
redshifts). We find 95% of the sources are well fitted by at 
least one of the galaxy templates. The templates are then 
used to determine the redshift at which the source would 
drop out of the sample as a consequence of both cosmologi- 
cal dimming and fc-correction, and this is used as the Vmax 
if it is found to be lower than the Vmax determined from the 
radio data alone. 

We conduct the V/V m axtest with and without those 
sources with poor fits to the multi-colour data. The inclu- 
sion of the sources with poorer fits leads to a decrease in the 
V/V max estimate by ~ 0.01 at L 3 2 5 mhz < 10 25 W Hz -1 sr -1 , 
and ~ 0.02 at higher luminosities, thus we do not deem it 
to be an important uncertainty in our analysis. As well as 
the statistical errors, due to the size of the sample which 
decrease as (12N)~^, we investigate three other sources 
of error. First, the redshift estimate for each radio source 
has an associated uncertainty which is principally depen- 
dant on the magnitude errors and hence is a function of 
redshift. We use Monte-Carlo methods to assess the size of 
this source of error. We derive 1000 samples drawn from the 
true sample based on Gaussian deviations of the reported 
redshift error. Second, we re-derive the V/V ma x statistic for 
the sample after correcting for the systematic photometric 
underestimate seen in Fig. Finally - as we discuss in §2 - 
we investigate the effect of contamination of mis-identified 
sources on the V/Vmax statistic by altering the search ra- 
dius of FIRST-SDSS DR1 matched sources. We re-derive 
the V/Vmax statistic for search radii of 2, 3 and 4 arcsec sep- 
arations. 
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Figure 3. The V/Vmax statistic for low-luminosity (-L325MHZ < 
10 27 W Hz" 1 sr" 1 ) radio sources. Crosses denote the mean 
V/ Vmax statistic at a given redshift provided in Table The 
shaded regions are the 68%, 95% and 99% (dark to light shad- 
ing respectively) confidence intervals based on Gaussian devia- 
tions about the reported redshift. The open triangles are the 
V/Vmax statistic derived after correcting for the systematic un- 
derestimate in the photometric redshifts. 



5 RESULTS 

In Fig. 01 and Table we show the mean V/Vmax statistic in 
8 bins of radio luminosity in the range 10 23 < L325MHZ < 
10 27 W Hz- 1 sr" 1 . For each bin we provide the mean redshift 
(z) . The Ict error bars are the statistical errors derived from 
the sample size. As we discussed above the V/V m ax statistic 
is affected by the errors in the estimates of the redshift. 
Overplotted in Fig. 01 are the results of the Monte Carlo 
simulations at 68%, 95% and 99% confidence (dark to light 
shading respectively) based on Gaussian deviations about 
the reported redshift. We also re-derive the V/Vmax statistic 
for the sample after correcting for the systematic underesti- 
mate seen in Fig. 01 These adjusted estimates are provided 
in the last column of Tableland are shown as open triangles 
in Fig. El 

The 5% of objects that cannot be satisfactorily fit with 
the galaxy templates change the V/V m ax estimate by 0.02 at 
most. An inspection of the confidence limits in Fig. 01 sug- 
gests that the errors that arise from the uncertainty in the 
redshift are less significant than the statistical errors. Fur- 
ther, the systematic photometric error leads to a change in 
<V/Vmax > of ~ 0.01. We note the latter two sources of er- 
ror have the combined effect of moving the V/V m ax statistic 
towards V/V m ax = 0.5 for the L325MHZ < 10 25 W Hz -1 sr -1 
sample and away from it for the higher luminosity sample. 

The only evidence (> 2cr) for a strong positive evolu- 
tion comes from the powerful radio sources with L325MHZ > 
10 25 W Hz" 1 sr" 1 . It has long since been established that 
these powerful, predominantly FRII sources, exhibit positive 
evolution (i.e. V/V m ax > 0.5). However, the lower-luminosity 
population does not follow this trend. In contrast to the 
strong evolution indicated by the FRII sources, the lower 
luminosity sample shows little evidence for evolution up to 
z ~ 0.6. This dichotomy is illustrated further in Fig. |1] and 
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Table 2. The banded V/Vmax statistic as a function of redshift 
split at L 32 5MHz = 10 25 W Hz" 1 sr" 1 . 
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Tabled where we show the banded V/V max statistic, which 
enables us to see the evolution in the V/ V ma x statistic as a 
function of redshift (see Avni & Bahcall 1980) for two bins 
in radio luminosity [L325MHZ > 10 25 W Hz -1 sr -1 (upper) 
and L325MHZ < 10 25 W Hz -1 sr -1 (lower)]. This result re- 
mains unchanged if we alter the search radius of the optically 
matched FIRST and SDSS DR1 radio sources. 

Table |H] shows the same banded V/V max statistic in 
the L325MHZ < 10 25 W Hz" 1 sr" 1 bin as in Table H but 
for search radii of 4, 3 and 2 arcsec. The banded V/Vmax 
statistic remains consistent in each of the separations but 
the relative levels of completeness fall to 80%, 68% and 57% 
for 4, 3 and 2 arcsec separations respectively. However, due 
to the consistency in V/V ma x the completeness seems to be 
uniform in all 2 bins. 

The evidence for the lack evolution in the low- 
luminosity sources is broadly consistent with a radio sur- 
vey by Waddington et al. (2001) which uses 72 sources, 42 
of which have accurate spectroscopic redshifts. Examination 
of their Fig. 10 suggests no significant evolution of sources 
at < 2 < 1 which is consistent with our lowest luminosity 
bin of P 325 mhz ~ 10 23 W Hz -1 sr" 1 shown in Fig. El 

Jackson & Wall (1999) found that for 26 sources at 
the lower luminosity, <l / /V max > = 0.314±0.057. If we cal- 
culate the mean V/V max statistic for our low-luminosity 
and high-luminosity sources over all redshifts we find that 
the sample with 10 23 < L325MHZ < 10 25 W Hz" 1 sr" 1 
yields <V/V max > = 0.47±0.01 for 686 sources. The sam- 
ple L325MHZ > 10 25 W Hz" 1 sr" 1 has a <V/V max > = 0.58 
± 0.02. This is obviously calculated over a large range in 
redshift which may not be informative regarding the evolu- 
tion. Indeed two samples with different < V/Vmax > do not 
necessarily have different evolutions ( Anvi & Bahcall 1980) . 
However it does indicate that the low-luminosity sources are 
consistent with little or no evolution and that the Jackson & 
Wall (1999) values are ruled out at > 5<r. A direct compari- 
son between such studies and the one presented here is made 
difficult because of the small sample size so that the samples 
may be drawn from completely different distributions. 

We also note that our results are consistent with Snellen 
& Best (2001) as their two sources are at L325MHZ ~ 
10 25 W Hz" 1 sr -1 and as such are at the cross-over point 
between the lower and higher luminosity sources in our anal- 
ysis. 



Figure 4. The banded V/V max statistic as a function of red- 
shift for high-luminosity sources (L325MHZ > 10 25 W Hz -1 sr - 1 ; 
upper half) and low-luminosity radio sources (L 3 25MHz < 
10 25 W Hz -1 sr -1 ; lower half). In the plot V e is volume enclosed 
at the redshift of the source, V a is the volume which corresponds 
to the maximum redshift which the source can be seen and Vo 
is the volume which corresponds to zq. The shaded regions have 
the same meanings as Fig. 151 



6 SUMMARY AND DISCUSSION 

We have conducted the first large scale analysis of the evo- 
lution in the comoving space density of radio sources at 
the faint end of the radio luminosity function. We match 
steep spectrum sources from the WENSS, NVSS and FIRST 
radio catalogues to the SDSS DR1 photometric catalogue 
which results in 1014 sources with redshift information ac- 
curate to Az ~ 0.03 at 2 < 0.5. Of these sources 686 
sources have L325MHZ < 10 25 W Hz -1 sr -1 and 328 have 
-k325MHz > 10 25 W Hz -1 sr -1 . For the powerful sources, as- 
sumed to be predominantly FRII radio sources, we are in 
good agreement with previous work based on spectroscopi- 
cally complete radio samples (e.g. Dunlop & Peacock 1990; 
Willott et al. 2001). However, for the lower luminosity pop- 
ulation, we find that the comoving space density remains 
approximately constant with increasing redshift, and thus 
these objects do not follow the trend of their more powerful 
counterparts. 

This is in rough agreement with what one would ex- 
pect by combining the information from both the radio lu- 
minosity function derived from a brighter, complete flux- 
density limited sample (e.g. Willott et al. 2003) and the 
low- frequency source counts. We rule out a strong neg- 
ative evolution, as suggested by Jackson & Wall (1999) 
at > 5a. We also find no evidence for a significant in- 
crease in the comoving space density of the sources with 
^325MHz < 10 25 W Hz -1 sr -1 . However, it is difficult to 
assess the significance of any possible dichotomy in the evo- 
lution of the low-luminosity and the high-luminosity sources 
with our sample alone, as we probe completely different 
ranges in redshift, due to redshift - luminosity degeneracy 
inherent to flux-density limited samples. In comparison to 
the combined 3CRR, 6CE and 7CRS samples of Willott et 
al. (2001) which covers a combined area of ~ 4.3 sr with 
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Table 3. The banded V/V max statistic for the low luminosity 
sources (L325MHZ < 10 25 W Hz -1 sr — 1 ) as derived in Table 2 
but for FIRST-SDSS DR1 matching radii of 4, 3 and 2 arcsec. 
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0.35 


0.09 


0.25 


0.10 


0.23 


0.13 


0.8 
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virtually complete redshift information, and has many pow- 
erful radio sources at 2 < 0.5, our analysis shows that the 
lower-luminosity population is unlikely to have evolved in 
such a dramatic way. Whether this is due to the evolution 
in comoving space density correlating with radio luminosity 
or with FRI/FRII morphology is impossible to determine 
with the resolution of the FIRST survey used in this analy- 
sis. Access to further high-resolution radio data would fully 
resolve this issue. 
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